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Abstract 



that is tuned to the high-p t nuclear modification factor measured in 
central Au-Au collisions at ^s NN = 200 GeV. Once tuned, the model 
is consistently applied to all the high-p t observables at 200 GeV: the 
r S^ . centrality evolution of the nuclear modification factor, the suppression 

of the away-side jet-like correlations, and the azimuthal anisotropies 
for these observables. Predictions for the leading-particle suppression 
at nucleon-nucleon centre-of-mass energies of 62.4 and 5500 GeV are 
presented. The limits of the eikonal approximation in the BDMPS 
approach, when applied to finite-energy partons, are discussed. 
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1 Introduction 

High-momentum leading-particle suppression in nucleus-nucleus (AA) with 
respect to proton-proton collisions is regarded as one of the major discoveries 
at the Relativistic Heavy Ion Collider (RHIC), Brookhaven. In Au-Au col- 
lisions at centre-of-mass energy */s NN = 200 GeV per nucleon-nucleon (NN) 
pair, the two experiments with high transverse momentum, p t , capabilities, 
PHENIX and STAR, have measured: 

• the suppression of single particles at high p t ( > 4 GeV) and central 
pseudorapidity (|?7| ^ 1), quantified via the nuclear modification factor 

p ( \ - l w d 2 A AA /dp t dr/ 

RAA(Pt) = JW~\ X A2AT I A A ' X 

(A coll ) ccntralityclass d 2 N pp /dp t dr] 

which would be equal to unity if the AA collision was a mere superposi- 
tion of N co \\ independent NN collisions (A^on scaling); instead, at high 
Pt -Raa is found to decrease from peripheral to central events, down 
to ~ 0.2 in head-on collisions [1,2]; the suppression is the same for 
charged hadrons and neutral pions for p t > 5 GeV; 

• the disappearance, in central collisions, of jet-like correlations in the 
azimuthally-opposite side of a high-p t leading particle [3] ; 

• the absence of such effects in d-Au collisions at the same energy [4, 5] . 

These observations can be naturally explained in terms of attenuation 
(quenching) of energetic partons produced in initial hard scattering processes, 
as a consequence of the interaction with the dense QCD medium expected to 
be formed in high-energy heavy-ion collisions. Several theoretical works exist 
on the subject [6-13]. Most of them implement the idea of parton energy 
loss due to medium-induced gluon radiation. 

In our Monte Carlo program PQM (Parton Quenching Model) we com- 
bine a recent calculation of parton energy loss [13] and a realistic description 
of the collision geometry, which was proven to play an important role [14]. 
Our approach allows to study and compare to RHIC data the transverse 
momentum and centrality dependence of single-hadron and di-hadron corre- 
lation suppressions, as well as the 'energy-loss induced' azimuthal anisotropy 
of particle production in non-central collisions. The model has one single 
parameter that sets the scale of the energy loss. Once the parameter is fixed 
on the basis of the data at a/s^ = 200 GeV, we scale it to different energies 
assuming its proportionality to the expected volume- density of gluons, as ar- 
gued in Ref. [15]. We then apply the same approach to calculate the nuclear 
modification factors at intermediate RHIC energy, ^/s NN = 62.4 GeV, and at 
LHC energy, ^/s NN = 5.5 TeV. Since we do not include so-called initial-state 



effects, such as nuclear modification of the parton distribution functions and 
parton intrinsic transverse-momentum broadening, we restrict our study to 
the high-pt region, above 4-5 GeV at RHIC energies and above 10 GeV at 
LHC energy, where these effects are expected to be small (less than 10% on 
Raa) [16,17]. 

2 Parton energy loss and collision geometry 

For the calculation of in-medium parton energy loss, we use the quenching 
weights in the multiple soft scattering approximation, which were derived 
in Ref. [13] in the framework of the 'BDMPS' (Baier-Dokshitzer-Mueller- 
Peigne-Schiff) formalism [7]. 

In a simplified picture, an energetic parton produced in a hard collision 
undergoes, along its path in the dense medium, multiple scatterings in a 
Brownian-like motion with mean free path A, which decreases as the medium 
density increases. In this multiple scattering process, the gluons in the parton 
wave function pick up transverse momentum k t with respect to its direction 
and they may eventually decohere and be radiated. 

The scale of the energy loss is set by the characteristic energy of the 
radiated gluons 

u c = qL 2 /2, (2) 

which depends on the in-medium path length L of the parton and on the 
BDMPS transport coefficient of the medium, q. The transport coefficient 
is defined as the average medium-induced transverse momentum squared 
transferred to the parton per unit path length, q = (k 2 ) mcdium /\ [13]. For a 
static medium it is time-independent. 

Differently from the original BDMPS calculation [7], in Ref. [13] the trans- 
verse momentum k t of a radiated gluon is kinematically bound to be smaller 
than its energy uj. The constraint k t < uj is imposed via the dimensionless 
quantity 

*=g = ^'. (3) 

which relates the scale of uj 2 , given by the square of the characteristic energy 
uj 2 , to that of k 2 , given by qL, as easily seen from the definition of q. The 
BDMPS case corresponds to R — > oo and it can be recovered by considering 
an infinitely-extended medium (L — > oo for fixed, finite, uj c ) [13]. 

The two parameters uj c and R determine the energy distribution of radi- 
ated gluons, uj dl/duj. While uj c sets the scale of the distribution, R controls 
its shape in the region < uj <C uj c , where the kinematic bound k t < uj is 



relevant. In the limit R — ► oo the distribution is of the form [7]: 

dl ^2a s C R f ^ for 



lim w r_ ~ -^^ ^ v 2- - ( 4 ) 

where Cr is the QCD coupling factor (Casimir factor) between the considered 
hard parton and the gluons in the medium; it is CV = 4/3 if the parton is a 
quark and Ca = 3 if the parton is a gluon. 

In the eikonal limit of very large parton initial energy E (E ^> u) c ), 
the integral of the radiated-gluon energy distribution estimates the average 
energy loss of the parton: 

f°° dl 

(AE) R = lim / w— duj oc a s C R uj c oc a s C R qL 2 . (5) 

R— >oo Jq O.0J 

Note that, due to the steep fall-off at large uj in Eq. (J3J), the integral is dom- 
inated by the region uj < uj c . The average energy loss (AE) is: proportional 
to a s C R and, thus, larger by a factor 9/4 = 2.25 for gluons than for quarks; 
proportional to the transport coefficient of the medium; proportional to L 2 ; 
independent of the parton initial energy E. It is a general feature of all 
parton energy loss calculations [7,9-11,13] that the radiated-gluon energy 
distribution uj dl/du does not depend on E. Depending on how the kine- 
matic bounds are taken into account, the resulting AE is .E- independent 
(BDMPS) [7] or depends logarithmically on E [9-11]. However, there is al- 
ways a stronger intrinsic dependence of the radiated energy on the initial 
energy, determined by the fact that the former cannot be larger than the 
latter, AE < E. Within the above simplified derivation which agrees with 
the main features of the BDMPS formalism, this kinematic constraint could 
be partially included by truncating the gluon energy distribution uj dl /du at 
the parton energy E. This would give (AE) oc a s C R ^fu~T c ^/min(u; c , E). For 
E < co c , we have (AE) oc a s C R \fqyE L: the kinematic constraint turns the 
L-dependence from quadratic to linear 1 and introduces a v^-dependence. 
Note that this procedure implements the constraint uj < E for the emission 
of one gluon, but it does not prevent from having AE = uji+uj 2 + ... > E in 
a multiple-gluon emission. A full theoretical treatment of the finite parton 
energy case in the BDMPS framework is at present not available. As we will 
discuss in Section El this introduces significant uncertainties in our results. 

The probability P(AE) that a hard parton radiates the energy AE due 
to scattering in spatially-extended QCD matter is known as the quenching 

1 DifFercnt approaches [9-11] emphasize the quadratic dependence of energy loss on the 
size of the medium down to rather small parton energies. 



weight [18]. In Ref. [13] the weights are calculated on the basis of the BDMPS 
formalism for quarks and gluons as a function of the two parameters u c and 
R and they are given as: 

P(AE; R, u c ) = Po (R) S(AE) + p(AE; R, u c ) . (6) 

The discrete weight po = Po(R) is the probability to have no medium-induced 
gluon radiation and the continuous weight p(AE) = p(AE; R, uo c ) is the 
probability to radiate an energy AE, if at least one gluon is radiated. In 
this work we use the quenching weights calculated in Ref. [13] with a fixed 
value of the strong coupling a s = 1/3. Note that, since (AE) oc a s q, the 
dependence on the value of a s can be largely absorbed in a rescaling of q. 

It has been shown [13] that a simple scaling law exists, which translates 
the radiated-gluon energy distribution for an expanding medium with a time- 
decreasing q(t) into an equivalent distribution for a static medium, with a 
time-averaged (q) = constant, using 



ti) = j2j (*-&) ?(£)<*, (7) 



where £o ~ 10" 1 fm <C L is the formation time of the expanding system. 

Due to the fact that q and L are two more intuitively and physically 
meaningful parameters, in all the previous applications [13,21,22] the natural 
(R, u; c )-dependence of the quenching weights was 'translated' into a (q, In- 
dependence, via Eqs. © and ©• The standard approach was to fix a value 
for the transport coefficient, the same for all produced partons, and then ei- 
ther use a constant length [13] or calculate a different length for each parton 
according to a description of the collision geometry [21,22]. However, this 
approach is not optimal, because (a) there is no unique and exact defini- 
tion of the in-medium path length when a realistic nuclear density profile is 
considered, as pointed out in Ref. [21], and (b) the medium density is not 
constant over the whole nucleus-nucleus overlap region but rather decreasing 
from the centre to the periphery. 

In order to overcome these limitations, we adopt a new approach in which 
the two parameters u c and R that determine the quenching weights are com- 
puted on a parton-by-parton basis, taking into account both the path length 
and the density profile of the matter traversed by the parton. 

Starting from Eq. (J2J) and using Eq. with a space-point dependent 
transport coefficient g(£) and a path-length-averaged (q), we define the ef- 
fective quantity 

1 f°° 

w c [effective = -(q) L 2 = / £q(£) d£ , (8) 



which on the r.h.s. does not explicitly depend on L. For a step-function 
'density' distribution g(£) = %9{L — £), Eq. (jHJ) coincides with Eq. (J2J). 
Similarly, we define 

/•oo 

(q) L (effective = / q(0 d£ (9) 

Jo 

and 

j-, I ^ V^c I effective ) / 1 n \ 

-ft I effective = / -\ r l ' \^J 

\q) L* | effective 

Using the definitions in Eqs. (|H)l- (|l()|l we incorporate the collision geometry 
in the calculation of parton energy loss via the 'local' transport coefficient 

3 Leading-particle suppression procedure 

Within the perturbative QCD collinear factorization framework, the expres- 
sion for the production of high-p t hadrons at central rapidity, y = 0, in pp 
collisions (no energy loss) reads: 



2„h 



dV 



dp t dy 



^ / dx a dx b dZj fa(x a ) fo(xb) 



^frab-tjX 



D h/j(Zj] 



z 2 

Vi=0 o 



(11) 

where f a ^) is the parton distribution function for a parton of type a(b) car- 
rying the momentum fraction x a ( b ), <r ab ~* jX are the partonic hard-scattering 
cross sections and Dh/j(zj) is the fragmentation function, i.e. the probability 
distribution for the parton j to fragment into a hadron h with transverse 
momentum p t = Zjp t j. To simplify the notation, we have dropped the de- 
pendence of a ab ^i x on y/s and of f a (b), a ab ~* : > x and Dh/j on the square of 
the scale (momentum transfer) Q 2 of the hard scattering, usually Q 2 ~ p 2 
Medium-induced parton energy loss is included by modifying Eq. (fTTj) to: 
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(12) 

where P(AEj]Rj,u c j) is the energy-loss probability distribution, Eq. (@), 
for the parton j (we will explain in the following how the input parameters 
R and u c for a given parton are calculated). 

In PQM we obtain the leading-particle suppression in nucleus-nucleus col- 
lisions by calculating the transverse momentum distributions in Eqs. (jll|l 



and (|T2"j) in a Monte Carlo approach. The 'event loop' that we iterate is the 
following: 

1. generation of a parton, quark or gluon, with p t > 3 GeV, using the 
PYTHIA event generator [23] in pp mode with CTEQ4L parton dis- 
tribution functions [24]; the p t -dependence of the quarks-to-gluons ratio 
is taken from PYTHIA; 

2. determination of the two input parameters, uj c and R, for the calcula- 
tion of the quenching weights, i.e. the energy-loss probability distribu- 
tion P(AE); 

3. sampling of an energy loss AE according to P(AE) and definition of 
the new parton transverse momentum, p t — AE; 

4. (independent) fragmentation of the parton to a hadron using the leading- 
order Kniehl- Kramer-Potter (KKP) fragmentation functions [25]. 

Steps 2 and 3 are explained in detail in the following paragraphs. Quenched 
and unquenched p t distributions are obtained including or excluding the third 
step of the chain. The nuclear modification factor i?AA(pt) is given by their 
ratio. Our hadrons pt distribution without energy loss at ^/s = 200 GeV 
agrees in shape with that measured for neutral pions in pp collisions by 
PHENIX [26]. 

Determination of u c and R 

We define the collision geometry in the (x, y) plane transverse to the beam 
direction z, in which the centres of two nuclei A and B colliding with an im- 
pact parameter 6 have coordinates (—6/2, 0) and (6/2, 0), respectively. Using 
the Glauber model [19] to describe the geometry of the collision, we as- 
sume (a) the distribution of parton production points in the transverse plane 
and (b) the transverse density of the medium both to be proportional to 
the 6-dependent product T^T B (x, y\ 6) = Ta(x, y) x T B (x, y) of the thickness 
functions of the two nuclei. The nuclear thickness function is defined as the z- 
integrated Wood-Saxon nuclear density profile: T[(x,y) = J dz p^ s (x,y,z). 
The parameters of the Wood-Saxon profile for different nuclei are tabulated 
from data [20]. Note that TaTb(x, y;b) estimates the transverse density of 
binary NN collisions, p C o\i(x,y; 6), modulo the inelastic NN cross section. 

We consider only partons produced at central rapidity and assume that 
they propagate in the transverse plane (E ~ p ~ p t ). For a parton with 
production point (x ,2/o) and azimuthal propagation direction (cos0o,sin0 o ) 
(</>o is sampled uniformly), we define the 'local' transport coefficient along 
the path of the parton inside the overlap region of the nuclei as: 

<?(£;6) = k x T A T B (x + £ cos (p ,y + £ sin O ; 6) , (13) 



where k is a free parameter (in fm) that sets the scale of the transport 
coefficient (in GeV 2 /fm). We compute the two integrals J and I\ (Eqs. 
and ©) 

/■oo 

I n = / emb)d£ n = 0, 1, (14) 



JO 

which determine the energy-loss probability distribution P(AE) using 2 (see 
Section EJ): 

u c = h and R=2I*/I . (15) 

Our approach allows a natural extension from central to peripheral nucleus- 
nucleus collisions: the idea is to fix the only free parameter, k, in order to 
describe the measured nuclear modification factor in central collisions and 
then use the impact parameter (6) dependence of the product Ta_Tb(x, y; b). 
We translate, by means of the Glauber model, the experimental definition 
of the centrality classes in terms of fractions of the geometrical cross section 
to a range in b and, within such range, we sample, for every loop of the 
chain reported at the beginning of this section, a value of b according to the 
6-dependence of the average number of binary collisions, d (N^/ /db. 

In order to give a synthetic and direct illustration of our results, we com- 
pute, for a given centrality class, the distributions of the two more customary 
variables L and q. To this purpose, for every parton we combine uo c and R, 
using Eqs. (J2j) and ®, to obtain an effective path length and an effective 
transport coefficient: 

L = R/u c = 2 Ji/Jo and q = 2u 2 J(LR) = / 7( 2/ i) • ( 16 ) 

We point out that the resulting definition of L is, as necessary, independent 
of k. Furthermore, it is the same one of us (A.D.) used in Ref. [21]. Note that 
q is proportional to k. In Fig. ^ we report for illustration the distributions 
of parton production points in the transverse plane and of in-medium path 
lengths, in central (0-10%), semi-central (20-30%) and peripheral (60-80%) 
Au-Au collisions. We will show the q distributions for different centralities 
in the next section (in Fig. |3J), after extracting the scale k from the data. 

Energy-loss sampling 

In the third step of the chain, we use the numerical routine provided in 
Ref. [13] for fixed a s = 1/3 to obtain the energy-loss probability distribution 
for given uj c , R and parton species (quark or gluon). According to this 
distribution, we sample an energy loss AE to be subtracted from the parton 

2 For simplicity, hereafter we drop the subscript "effective" for cu c and R. 
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Figure 1: Distributions of parton production points in the transverse plane 
(upper row) and in-medium path length (lower row) in central, semi-central 
and peripheral Au-Au collisions. The quantity (L) is the average of the path 
length distribution. 



transverse momentum. The quenching weights are calculated in the eikonal 
approximation, where the energy of the parton is infinite (E = p t = oo). 
Therefore, when the realistic case of finite-energy partons is considered, a 
significant part of the energy-loss probability distribution P(AE) lies above 
the parton energy E, in particular for large values of uj c and R, or equivalently, 
of q and L. The energy loss, under the constraints introduced by the finite 
parton energies, is sampled following two approaches: 

• Reweighted: truncate P(AE) at AE = E and renormalize it to unity by 
dividing out the factor J deP(e). The Monte Carlo implementation 
of this approach is: sample AE from the original P(AE); if AE > E, 
sample another AE; iterate until a AE < E is sampled. 

• Non-reweighted: truncate P(AE) at AE = E and add the (^-function 
5{AE—E) J E de P{e) to it. The integral of P is, in this way, maintained 
equal to one. The corresponding Monte Carlo implementation reads: 
sample an energy loss AE and set the new parton energy to zero if 
AE> E. 



The resulting energy loss is larger in the non-reweighted case, where par- 
tons are 'absorbed' by the medium with a probability j E deP(e). As we 
will see in the next section, the difference can be quite large for low pt and 
sufficiently-large transport coefficients. It is argued [13, 27] that the dif- 
ference in the values of the observables for the two approaches illustrates 
the theoretical uncertainties. Along the lines of what is done in a recent 
work [22] developed in parallel to the present study, we display our model re- 
sults as a band delimited by a solid line representing the non-reweighted case 
(larger quenching) and a dashed line representing the reweighted case (smaller 
quenching). For the time being, from the theory side both approaches are 
highly disputable, while the guidance given by the experimental results will 
be commented in the conclusions. 



4 Results 

Nuclear modification factor in Au— Au at ^/s NN = 200 GeV 

We start by presenting the results on RAA^Pt) in central Au-Au collisions 
at y/s NN = 200 GeV obtained using constant in-medium path length and 
transport coefficient (left-hand panel of Fig.EJ). The data on charged hadrons 
and neutral pions from PHENIX [1] and STAR [2] are reported with combined 
statistical and p t -dependent systematic errors shown by the bars on the data 
points and p t -independent normalization errors shown by the bars centred at 
-Raa = 1- The model results are shown by the lines: for all hadrons, with 
q = 1 GeV 2 /fm and L = 6 fm, (solid line) and for hadrons coming from 
quarks and from gluons, separately, with q = 0.75 GeV 2 /fm and L = 6 fm, 
(dashed and dot-dashed lines). In order to compare our results to those in 
Ref. [13], we use the same parameters and treat the finite-energy constraint 
in the non-reweighted case. The two lines obtained with q = 0.75 GeV 2 /fm 
and L = 6 fm agree with those reported in Fig. 20 of Ref. [13]. Since the 
high-p t hadron spectrum at RHIC energies is mainly coming from quarks 3 , 
which lose less energy than gluons, a larger q of ~ 1 GeV 2 /fm is necessary 
to match the measured -Raa, when a realistic quarks-to-gluons ratio is used. 
Before moving to the parton-by-parton approach of PQM outlined in the 
previous sections, it is very instructive to show the model results obtained 
using a constant transport coefficient and the Glauber-based path-length 

3 At y/s = 200 GeV, with CTEQ4L parton distribution functions [24], gluons dominate 
the parton p t distribution up to about 20 GeV. However, since quarks fragment harder 
than gluons, high-p t hadrons are mostly produced from quark fragmentation. Using KKP 
fragmentation functions [25] , we find that 75% of the hadrons with p t > 5 GeV come from 
quark fragmentation and 25% from gluon fragmentation. 
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Figure 2: i?AA(pt) for central Au-Au collisions at ^/s NN = 200 GeV. 
PHENIX [1] and STAR [2] data are reported with combined statistical and 
p t -dependent systematic errors (bars on the data points) and pt-independent 
systematic errors (bars at -Raa = !)• Model results for constant q and L (left- 
hand panel) and for constant q and Glauber-based L distribution (right-hand 
panel) are reported. In the right-hand plot and in all the following figures the 
shaded band is delimited by non-reweighted case (solid line) and reweighted 
case (dashed line). 



distribution for 0-10% central collisions (displayed in the bottom-left panel 
of Fig. EJ). We have to use q ~ 15 GeV 2 /fm to describe the data with the 
model band delimited by the reweighted and non-reweighted cases (right- 
hand panel of Fig. EJ). When going from a constant L — 6 fm to a realistic 
distribution, the transport coefficient has to be increased by more than one 
order of magnitude, because there are many partons with small path lengths 
of 2-3 fm that can be quenched only if the medium is very dense. It is 
interesting to note that Raa is clearly increasing with p t when a constant 
length is used, while it is flatter with the full distribution. This is due to the 
presence of a long tail in the L distribution, up to 12 fm: only high-energy 
partons can fully 'exploit' this tail, while low-energy ones are just completely 
stopped by the medium after a few fm, so that the 'effective' average length 
increases with the parton energy. We note that between the non-reweighted 
and reweighted approach to the parton finite-energy constraint there is a 
difference of about a factor 2 in the magnitude of Raa, but also a difference 
in the slope versus p t , which is slightly positive for non-reweighted and slightly 
negative for reweighted. 

Using a constant transport coefficient of 15 GeV 2 /fm and a realistic L 
distribution, the measured hadron suppression can be fairly well described for 
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Figure 3: i?AA(Pt) for central Au-Au collisions at ^ s nn = 200 GeV. The 
model band is obtained with a parton-by-parton calculation of uj c and R. 
The average transport coefficient is 14 GeV 2 /fm. 



p t ~ 5 GeV (at lower p t we do not apply the model, as initial-state effects and 
in-medium hadronization, that we do not include, might play an important 
role). Remarkably, our result agrees with that obtained in Ref. [22], where 
the same quenching weights and a simplified collision geometry with effective 
nuclei (cylindrical density profile instead of the Wood-Saxon we use) are 
coupled to a leading-order perturbative QCD calculation. The R\a band is 
found to have similar p t -dependence (rather flat) and width. Numerically, 
the extracted value of q is ~ 10 GeV 2 /fm in Ref. [22], smaller than our 
15 GeV 2 /fm. However, this is not an inconsistency, since the value of a s 
used in the calculation of the quenching weights is 1/2 in Ref. [22] and 1/3 
here, and the scale of the energy loss is set by the product a s q (see Eq. (J3J)). 
In order to address the centrality dependence of the high-p t suppression, 
we move to the parton-by-parton approach. For central collisions, the re- 
sult obtained with the scale parameter k = 5 x 10 6 fm, corresponding to 
(q) ~ 14 GeV 2 /fm, is shown in Fig. El The model band is very similar to 
that reported in the right-hand panel of Fig. El for q = 15 GeV 2 /fm and 
the L distribution. We now vary the centrality, keeping always the same 
scale k. Figure 0] shows the distributions of q, calculated from Eq. (|TB|l . for 
different centrality bins. The q variation within a given bin reflects the differ- 
ent parton production points, hence different medium densities encountered. 
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Figure 4: Distributions of q, calculated from Eq. (fTB|) . for different central- 
izes; the k parameter is fixed to the value that allows to describe -Raa for 
the most central collisions. 




Figure 5: -Raa(p0 for different centralities. Data are PHENIX charged 
hadrons (closed squares) and tt° (open squares) [1] and STAR charged 
hadrons (stars) [2]. 
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Figure 6: Average Raa in the range 4.5 < p t < 10 GeV [1,2] (left-hand panel) 
and /aa, defined in the text, for the away-side jet [3] (right-hand panel) as 
a function of collision centrality, expressed by the number of participants, 
-Wpart- F° r -Raa, the error bars are combined statistical and p t -dependent 
systematic errors and the bands centred at Raa — 1 are the p t -independent 
normalization errors for PHENIX (dashed) and STAR (dot-dashed). For 
/aa, the statistical (bars) and systematic (ticks) errors are shown. 



The rightmost (highest) value refers to partons originating in the centre of 
the collision system. The model nuclear modification factors, compared to 
PHENIX [1] and STAR [2] data ranging in centrality from 0-5% to 80-92%, 
are reported in Fig. 0D We note that the theoretical uncertainty band is 
narrower for semi-central and peripheral collisions, where, due to the smaller 
size and density of the medium, the probability to have AE > E in the 
quenching weights becomes marginal. 

Our results follow the decrease of the measured Raa with increasing cen- 
trality. This is more conveniently visualized in the left-hand panel of Fig. 
where we show the average -Raa in the range 4.5 < p t < 10 GeV plotted as 
a function of the number of participant nucleons, iVp ar t, obtained from the 
Glauber model. Data are taken from Refs. [1,2]. 



Back-to-back correlations 

By generating pairs of back-to-back partons, we can study the centrality 
dependence of the disappearance of the away-side jet. This effect is usually 
quantified using the correlation strength [28] 



D 
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dp 



t.2 



d 3 rr hlh2 
dA0- AA 



/dp t ,idpt,2dA0 



A0>A?i n 



d^AA/dPt, 



(17) 
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for an associated hadron h 2 with transverse momentum p t}2 in the opposite 
azimuthal direction of a trigger hadron hi with transverse momentum p tl . 
The STAR data [3] are for trigger particles with 4 < p t ,i < 6 GeV and 
associated particles with p t<2 > pT m = 2 GeV and p t ,2 < Pt,i, with A0 = 
|0i — 02 1 > A0 min = 130°. The correlation strength is then corrected for 
combinatorial background and azimuthal anisotropy of particle production in 
non-central collisions [3]. The correlation strength in nucleus-nucleus relative 
to pp collisions defines the suppression factor: 

/aa = ^. (18) 

We generate pairs of partons with the same initial p t and separated in azimuth 
by A0 = 180°. Then, we calculate uj c and R for each parton and apply energy 
loss and fragmentation. We count as trigger particle every hadron h\ with 
4 < p t i < 6 GeV and as associated away-side particle the other hadron h 2 
of the pair, if its transverse momentum is in the range 2 GeV < pt,2 < Pt,i- 
We define: 



'AA 



lyassociatcd \ / / ?yassociated 



yytrigger / / V yytrigger 

' with energy loss I \ / w/o energy loss 



(19) 



The right-hand panel of Fig. |U1 shows our result for Jaa versus N pSuTt , com- 
pared to STAR measurements in Au-Au collisions at ^/s NN = 200 GeV, with 
statistical (bars) and systematic (ticks) errors, from Ref. [3]. The magnitude 
and centrality dependence of the suppression are described without changing 
the scale parameter value we extracted from -Raa in central collisions. 

Azimuthally-differential observables 

For non-central collisions, the nucleus-nucleus overlap profile is asymmetric 
with respect to the event-plane direction, defined by the line that contains 
the centres of the two colliding nuclei and the impact parameter vector b, in 
the transverse plane. The asymmetry is visible in the upper row of Fig. ^ 
where the event-plane direction is parallel to the x axis. Consequently, the 
in-medium path length is, on average, larger for partons propagating in the 
out-of-plane direction (perpendicular to the event plane) than for partons 
propagating in the in-plane direction (parallel to the event plane). 

Due to parton energy loss, the asymmetry in the medium geometry should 
be reflected in the azimuthal distribution dA/d0 of high-p t hadrons with 
respect to the event plane, = 0°. We quantify this effect by calculating: 

• the azimuthal anisotropy, as given by the second Fourier coefficient of 
the dN/d(j) distribution, v 2 [29] ; we obtain the value of v 2 for hadrons 
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Figure 7: Transverse momentum and centrality dependence of the az- 
imuthal anisotropy 1)2, compared to measurements for charged hadrons from 
PHENIX [30] and STAR [31,33] (see text). Only the statistical errors are 
plotted. Note that, opposed to all other figures, here the model result for the 
non-reweighted case (solid) is the upper limit of the band and that for the 
reweighted case (dashed) is the lower one. 



in a given p t range by fitting their azimuthal distribution to the form 

a ■ (1 + 2 v 2 cos 20); 

fi^(pt), the nuclear modification factor for hadrons in an azimuthal 



cone of 45° centred at the angle 
we use cf>o = 0° (in-plane), O = 
(intermediate); 

r0O 



bo with respect to the event plane; 
90° (out-of-plane) and <p = 45° 



(away side), the nucleus-nucleus away-side correlation strength 
relative to pp, in the three azimuthal regions defined for R^a- 

The scale parameter k is again kept to the value that allows to match the 
measured .Raa m central collisions at ^s NN = 200 GeV. 

Figure U\ (left-hand panel) shows the model results for v 2 as a function 



of the transverse momentum, for central (0-10%, N ] 
central (20-60%, JV, 



part 



part ~ 320) and non- 
100) Au-Au collisions, compared to non-central 
experimental measurements on charged hadrons obtained by PHENIX [30] 
and STAR (preliminary) [31] using three different methods: reaction plane 
reconstruction (t^jRP}), 2-particle correlations (1*2 {2}) [32] and 4-particle 
correlations (i>2 {4}) [32]. In the right-hand panel of the same figure, the v 2 
centrality dependence from the model is compared to charged hadrons data 
from PHENIX [30], at p t w 4.5 GeV, and from STAR (preliminary) [33], at 
p t ps 6 GeV. 

The measured azimuthal anisotropy at intermediate transverse momenta 
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Figure 8: Azimuthal variation of the nuclear modification factor of the away- 
side correlations. In the left-hand panel, ratios Ra°a/Raa (out-of-plane/in- 
plane) and Raa/R-aa (intermediate/in-plane), averaged over the range 4.5 < 
p t < 10 GeV. In the right-hand panel, the same ratios for /aa with trigger 
conditions as in Fig. [HI Both observables are plotted as a function of N pauTt . 



of 4-6 GeV is systematically larger than that generated by parton energy loss 
in our model, indicating the presence of non-negligible collective flow effects 
in this momentum range. However, the preliminary STAR measurements 
at higher p t , shown in the left-hand panel of Fig. suggest that v 2 might 
go down to values compatible with those expected from parton energy loss 
in an azimuthally-asymmetric medium. High-p t data with larger statistics 
from the recent RHIC Run-4 will allow to clarify this point. We note that 
our maximum v 2 of 0.05-0.10, for A^ part « 100, is similar to that obtained in 
other parton energy loss [28] or absorption [14] calculations. 

The azimuthal variation of the nuclear modification factor and of the 
away-side correlations is illustrated in Fig. |S1 For R^a (left-hand panel), 
we show the two ratios R%°a/Raa (out-of-plane/in-plane) and R\\/ R\a 
(intermediate/in-plane), averaged over the range 4.5 < p t < 10 GeV, as 
a function of collision centrality (N paTt ). As for v 2 , the asymmetry is maxi- 
mum at A^part ~ 100, where the model gives for Raa a ratio out-of-plane/in- 
plane of ~ 0.75. Similarly, for the away-side correlation Taa (right-hand 
panel), we show the two ratios I^a/Iaa (out-of-plane/in-plane) and I^a/^aa 
(intermediate/in-plane). The conditions on the near-side trigger and the as- 
sociated away-side particles are the same as for Fig. H3 At A part « 100-150 
the model predicts an away-side correlation strength of about 30% lower for 
the out-of-plane relative to the in-plane direction. Both effects are rather 
strong and their measurement at RHIC would be of great interest. 
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Figure 9: Model results for -Raa(p0 i n central and semi-peripheral Au-Au 
collisions at ^/s NN = 62.4 GeV. The preliminary ir° data (0-10% centrality 
class) from PHENIX [35] are also shown; the pp reference is the PHENIX 



PP 



7T 



X parameterization, the error bars on the data points are the 
combined statistical and ^-dependent systematic errors and the bar centred 
at Raa = 1 is the systematic error on the normalization. 



Nuclear modification factor at 



62.4 GeV 



The recent RHIC run with Au-Au collisions at -y/s NN = 62.4 GeV allows 
the measurement of the nuclear modification factor for charged hadrons and 
neutral pions up to transverse momenta of 7-8 GeV. We estimate the leading- 
particle suppression due to parton energy loss at this lower centre-of-mass 
energy by using the proportionality of the transport coefficient q to the ini- 
tial volume-density of gluons n gluons [15]. In the saturation model [34], for 



collisions of two nuclei with mass number A at energy 
estimated to scale as 



??.' 



gluons 



ocA 



0.383 



v 0.574 



such density is 



(20) 



0.5 x n A ™Au,200GeV- Applying this scaling to 



mi • • gluons 

This gives n Au _ Au624GeV 
the value of the k parameter, see Eq. ()13j) . found in our model for central 
collisions at 200 GeV, we obtain a transport coefficient distribution with 
mean value (q) ~ 7 GeV 2 /fm in central collisions at 62.4 GeV. 

We generate hard partons using PYTHIA at ^/s = 62 A GeV and use the 
procedure described in Sectional The results are shown in Fig. El along with 
preliminary data from PHENIX [35] for neutral pions up to p t ~ 7 GeV in 
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Figure 10: Raa{Pi) for different centrality classes in Pb-Pb collisions at 
y/s NN = 5.5 TeV and in Au-Au collisions at y/s NN = 200 GeV. For compar- 
ison, the prediction for LHC presented in Ref. [17] (Vitev-Gyulassy) is also 
shown. 



0-10% central collisions. For p t > 5 GeV, we find -Raa — 0.3, in accordance 
with the data, in central (0-10%) and ~ 0.7 in semi-peripheral (40-50%) 
collisions. These values are not much larger than those at ^s NN = 200 GeV. 
At smaller ^s NN , although the transport coefficient is reduced by a factor 
of 2, the increased softness of the parton transverse momentum distribution 
determines a stronger effect of energy loss on the nuclear modification factor. 
Prior to the release of the preliminary PHENIX data, several predictions 
were published [36-38] . While the magnitude and p t -dependence of Raa in 
Ref. [36] seem to agree with our result (although Raa is presented only up 
to p t = 6 GeV there), the predictions in Refs. [37,38] show a different trend 
with p t : Raa for Pt ~ 5 GeV is decreasing with increasing p t , down to values 
of about 0.2 at p t ~ 16 GeV. 

Extrapolation to the LHC 

To compute the expected nuclear modification factor in Pb-Pb collisions 
at the LHC we use PQM with the parton p t distribution extracted from 
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PYTHIA at \/s = 5.5 TeV. Scaling the k parameter according to Eq. (|20p. 
we have n^b,5.5TeV ^ 7 x ng™™ u ,200GeV> i.e. (q) ~ 100 GeV 2 /fm. 

We report in Fig. EH the expected transverse-momentum dependence of 
i?AA in the range 10 < p t < 90 GeV for different centralities (the results at 
y / s NN = 200 GeV are shown as well). In the most central collisions -Raa is 
of ~ 0.15, independent of p t . This value is about a factor of 2 smaller than 
that measured at ^fs^ = 200 GeV. Our result for the LHC is in agreement, 
both in the numerical value and in the ^-dependence, with that obtained 
in Ref. [22] using the same quenching weights and the same a s (q), while it 
is quite different from that calculated in Ref. [17] assuming an initial gluon 
rapidity density dAf gluons /d?/ in the range 2000-3500. For comparison, we 
have reported in the same figure the result of Ref. [17]: there, -Raa is predicted 
to rise significantly at large transverse momenta, from 0.1-0.2 at 20 GeV to 
0.4-0.6 at 90 GeV. We note that the difference between the two results is not 
likely to be due to the fact that we do not include nuclear ( ant i-) shadowing 
effects, since these are expected to determine a rather p t -independent increase 
of R AA of about 10% in the range 25 < p t < 100 GeV [16, 17]. 

5 Discussion 

High-energy partons from the surface 

The centrality dependence of leading-hadron suppression and back-to-back 
di-hadron correlations is well described by our model, in which the centrality 
evolution is purely given by collision geometry. This suggests that the high- 
opacity medium formed in Au-Au collisions at y/s NN = 200 GeV has initial 
size and density that decrease from central to peripheral events according to 
the overlap profile of the colliding nuclei, Ta(x,i/) x T^{x, y). At the centre 
of the medium the density is maximum and partons crossing this region are 
likely to be completely absorbed. Only partons produced in the vicinity 
of the surface and propagating outward can escape from the medium with 
sufficiently- high energy to fragment into hadrons with more than few GeV 
in p t . Such an 'emission from the surface' scenario was pictured also in a 
recent work [14], where the centrality dependence of -Raa and Jaa could 
be reproduced by a simple model of parton absorption whose only physical 
ingredient was a Glauber-based nucleus-nucleus overlap profile. 

The region from which partons escape from the medium is visualized 
by plotting the distribution of production points for partons that give a 
high-energy hadron (ph adron > 5 GeV). This distribution for central Au-Au 
collisions at 62.4 and 200 GeV and Pb-Pb collisions at 5.5 TeV is shown in 
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Figure 11: Distributions of parton production points in the transverse plane 
(upper row) and in-medium path length (lower row) for partons that escape 
the medium and produce hadrons with p t > 5 GeV in central Au-Au colli- 
sions at 62.4 and 200 GeV and in central Pb-Pb collisions at 5.5 TeV. The 
quantity (L) is the average of the path length distribution. These plots were 
obtained in the non-reweighted approach. 



Fig. ^2 along with the corresponding path length distribution. The 'thick- 
ness' of the escape region is of order 2-3 fm and it decreases as ^s NN increases 
from intermediate RHIC energy to LHC energy. 

It is interesting to try to apply a simple toy model: all partons with 
a path length L smaller than a maximum length L^£ escape from the 
medium, the others are absorbed. We define the path length probability dis- 
tribution V(£) as the probability distribution for a generic parton to have a 
path length £. The distributions in the lower row of Fig. ^ are examples of 
V{£) for different centrality classes in Au-Au. V{£) is normalized to unity, 
/ °° d£V(£) = 1, and, thus, the integral f Q d£V(£) gives the fraction of par- 
tons with path length smaller than L. Using the measured (or expected) 
Raa for given collision energy and centrality, 

_ T max , k 

j^escapc diV ^ = Raa _ At v / i ^ = 200 GeV, 
central (0-5%) to semi-peripheral collisions (40- 



-^escapc can De estimated as 



we find L™* pe w 



2.5 fm from 



in this wide central- 
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a: both absorbed 
b: one escapes 
c: both escape 




Figure 12: Correlation between path lengths for pairs of partons produced 
back-to-back in 0-10% central Au-Au collisions. 



ity range the simultaneous decrease of system density and volume results 
in energetic partons being emitted from a shell of constant thickness. For 
more peripheral collisions the system becomes very diluted and partons can 
escape from the whole volume (£™^ pe ~ 3.5-4 fm « system size). In central 
collisions at different energies, we find L™^ pe ps 3 fm at ■ s /s tiN = 62.4 GeV 
and L™* pe « 1.5 fm at ^/s^ = 5.5 TeV (LHC). 

Remarkably, this absorption toy model allows to reconcile the magnitude 
of single-particle and away-side correlation suppressions measured in central 
Au-Au at 200 GeV, as illustrated in Fig. ^J The figure shows the distri- 
bution L\ versus L 2 for pairs of partons (1 and 2) generated back-to-back. 
Parton pairs produced in the middle of the overlap profile populate the cen- 
tral part of the distribution (Li ~ L 2 ), while pairs produced closer to the 
surface are in the two tails (Li 3> L 2 or L\ <C L 2 ). We report the two 
lines L\ = L™£ and L 2 = £™^„, which divide the distribution in three 



T max ,, „ J T T max 

-^escape ana ^2 — -^escape' 

max 
escape '•"•' v " ^^ w^u ^j.^ ^^^^j.^^^, ^^j ±^± ^ i( ^ ^ -"^escape 

> L™^ only one of the two partons escape the medium, and (c) 
-£wL„ both partons escape. With the value L™^L„ = = 2.5 fm, 



parts: (a) for Lip > £™ca Pe both partons are absorbed, (b) for Li( 2 ) < L 



< 



and L 

IOr L/l t 2 ■ ^escape """■"■ t"™- «""" v -"""'r v - ..*«." ^^^ .^^ ^ escape 

extracted from the measured -Raa, the third part of the distribution (c) is 
empty: it never happens that both partons can escape, in agreement with 
the value compatible with zero measured by STAR for Jaa- 
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Figure 13: Average relative energy loss versus parton energy for quarks and 
gluons in central collisions at RHIC and LHC energies for the non-reweighted 
and reweighted cases. 



Energy-loss saturation 

The strong parton absorption suggests that we are in a saturation regime 
of the energy loss, AE/E — > 1, as almost all hard partons produced in 
the inner core are thermalized (AE/E = 1) before escaping the medium. 
Indeed, the average relative energy loss, (AE/E) (from the Monte Carlo), 
shown versus parton energy E in Fig. ^] for central collisions at ^s NN = 200 
and 5500 GeV, is almost saturating to unity for gluons (70-80%) and it is 
very large also for quarks (50-70%). Due to the fact that gluons are closer 
to energy-loss saturation than quarks, the ratio of gluon to quark (AE/E) is 
much smaller than the Casimir ratio Ca/Cf = 2.25 expected from Eq. (0). 
Furthermore, since absorption (and, hence, saturation) is more significant 
for small-i? partons, or, in other words, large--E partons can exploit larger 
energy losses, the genuine BDMPS AE/E oc 1/E is replaced by a rather 
^-independent effective AE/E. It is important to point out here that the pt- 
independent nuclear modification factor obtained in our model, in agreement 
with RHIC data above ~ 5 GeV, is a natural consequence of this saturation 
scenario. 

As the average relative energy loss is close to one, at least for the non- 
reweighted case, we are not very sensitive to the shape of the continuous 
part of the quenching weights in Eq. (JHJ), p(AE). Rather, the energy- loss 
probability is dominated by the discrete part, the probability to have no 
medium-induced radiation, p Q . In order to confirm this statement, we repeat 
the calculation with a modified PQM version: in the quenching procedure 
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Figure 14: Average -Raa in the range 4.5 < p t < 10 GeV (left-hand panel) 
and Jaa for the away-side jet (right-hand panel) as a function of the number 
of participants. The data points and errors are the same as in Fig. |S1 Here 
the model results obtained using either zero or maximum energy loss, as ex- 
plained in the text, (labelled 'only P(AE = 0)', lighter band) are shown 
together with those obtained with the standard procedure (labelled 'full 
P(AE)\ darker band). 



we consider the parton as absorbed whenever the sampled energy loss AE is 
larger than zero. That is, we have either no energy loss or maximum energy 
loss. Also in this case, we consider the two finite-energy constraint methods, 
non-reweighted and reweighted. We note that the 'survival' probability is po 
for the non-reweighted case and Po/ J deP(e) for the reweighted case. In 
Fig. EUwe report, as a function of the number of participants, -Raa and Iaa> 
calculated with this modified quenching procedure, labelled 'only P(AE = 
0)', and the same value of k we used for the standard procedure. For the 
most central collisions, down to N part — 150, the agreement with data is 
very good, whereas deviations are clearly visible in Raa when going to semi- 
peripheral and peripheral collisions, iVp art < 150. This confirms that, in 
central collisions, partons are either completely absorbed or coming from 
the surface, whereas in non-central collisions the shape of the energy-loss 
probability distribution plays a role in the description of the data. 



6 Conclusions 

Most of the present high-momentum observables have been studied using 
the Parton Quenching Model (PQM) in which hard partons are generated 
with PYTHIA [23], medium- modified with the quenching weights [13], and 
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hadronized independently via KKP fragmentation functions [25]. Using a 
Glauber approach with Wood-Saxon density profiles of the colliding nuclei, 
the chain takes into account the realistic spatial distribution of hard parton 
production points and the amount and density of matter traversed by each 
parton. 

The results show that, if parton production coordinates and realistic den- 
sity profiles are taken into account, the ensuing transport coefficient has to 
acquire very large values: (q) ~ 14 GeV 2 /fm in central Au-Au collisions 
at ^/Snn = 200 GeV. Note that we used a relatively small value of a s 
(1/3). Due to the scaling AE oc a s q, if larger values (e.g. 1/2) were used, 
the extracted transport coefficient would be smaller, but still quite large, 
(q) ~ 10 GeV 2 /fm [22]. In Ref. [22], it is pointed out that such q values do 
not necessarily imply unexpectedly large medium initial energy densities (a 
few hundreds GeV/fm 3 ), as one obtains in the hypothesis of an ideal plasma 
whose constituents interact perturbatively with the hard partons [15], but 
rather suggest that the medium might interact with the hard partons much 
stronger than perturbatively expected. Technically, the large extracted q 
values present a yet unsolved theoretical problem, since the eikonal approach 
used in the theory cannot be flawlessly extended to finite (low) parton ener- 
gies. In the situation, we presented two possibilities: 

1. the theoretical treatment is applied regardless of the obvious problem 
for low-energy partons, which are completely absorbed with a rather 
large probability (non-reweighted approach); 

2. a reweighting procedure is performed in order to prevent the complete 
parton absorption in the medium (reweighted approach). 

Our calculations contain only one free parameter that was adjusted to the 
measured nuclear modification factor in central collisions at y/s NN = 200 GeV 
at RHIC. The same parameter was then employed to extract the centrality 
dependence of the nuclear modification factor for hadrons and di-hadrons, 
and the azimuthal anisotropy parameter V2 ■ Some of these observables were 
simulated for ^s NN = 62.4 GeV and 5.5 TeV collisions, scaling only the 
expected medium densities. 

When comparing to experimental results, we observe a tendency for the 
outright application of the theory, i.e. non-reweighted approach, to fit the 
i?AA and Taa centrality dependence reasonably well and, in general, better 
than the reweighted case. The calculated value of V2 lies about 2 standard 
deviations below the experimental data at p t — 4-6 GeV, suggesting the 
presence of collective elliptic flow effects up to these p t values; upcoming 
measurements with higher statistics should be able to give a definitive an- 
swer whether elliptic flow still subsists at momenta larger than 7-8 GeV 
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or not. We have also shown some predictions for the azimuthal variation 
of leading-particle suppression and jet-like correlations. We note that as a 
consequence of the azimuthal anisotropy for high-p t particles simulated in 
the present approach the low-energy particles 'radiated' by quenched par- 
tons might contribute a v% of the 'opposite sign' in the low-p t region. Thus, 
the parton quenching at the LHC could produce an apparent decrease of 
the elliptic flow at low p t with respect to the predictions of hydrodynamic 
calculations. 

We observe that the reweighting 'simulates' a softer transport coefficient, 
i.e. it de facto allows for some partons to be emitted from the away side, 
while a simple geometrical toy model excludes the away-side partons once the 
model has been tuned on the measured -Raa- Furthermore, in the reweighted 
approach, Raa(Pi) is larger towards low transverse momenta, see e.g. Fig. El 
because the survival probability for a parton with energy E, p / J Q deP(e), 
increases when E decreases. This feature appears to be unphysical, or at 
least non-intuitive. Clearly, the full treatment of the difficulties encountered 
here should be tackled theoretically in a more complete way. 

The inspection of the production-point distribution for energetic partons 
escaping the medium and of the average relative energy loss suffered by 
quarks and gluons in central collisions at top RHIC energy depicts the dense 
medium in the nuclear overlap region as a black disk: either the partons are 
absorbed or they escape from a thin shell close to the surface. 

The present model, applied to the LHC, gives the interesting result that 
the i?AA value is essentially constant with p t , and very low, up to the highest 
parton energies. As shown in Fig.ES this prediction differs substantially from 
others obtained for the LHC [17]. Namely, in our model the black disk effect, 
which requires a large transport coefficient, extends the strong suppression 
up to very high transverse momenta. This scenario would amount to decrease 
the number of high-energy jets by almost an order of magnitude and it should 
be considered in the future planning of experimental studies. 
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